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SUMMARY 


Axial -field permanent magnet motors have been shown 
to be particularly suitable for electric vehicle drives, 
as evidenced by the support for such advanced motors in 
the Department of Energy's Electric and Hybrid Vehicle 
Program. There were two most important factors in the 
analytical modeling of this type of machine that had not 
been extensively studied, and which were the subject of 
this research. These are: 

(i) The modeling of anisotropic alnico-type permanent 
magnets . It is well-known that magnets with linear 
demagnetization characteristics can be simply modeled 
by an equivalent solenoid for the purposes of field 
computation. However, those with non-linear de- 
magnetization curves are very much more difficult to 
model, since a non-uniform magnetization distribution 
exists throughout the magnet volume. This research 
was aimed at developing a realistic analytical model 
for such materials, which predominantly include the 
alnico types. 

(ii) Modeling of the mechanical commutation in these motors . 
The high coercive force of permanent magnets allows 
large air-gaps to be tolerated, and this leads to the 
use of iron-less disc-type armatures and the axial-field 
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topology. The advantages aie the elimination of eddy-current 
and hysteresis power losses, and the gaining of high motor 
power densities. However, the consequent field in the 
magnetic neutral regions and the low inductances of the 
armature coils leads to greater short-circuit currents during 
commutation. The predominant power loss in this machine is 
therefore associated with the mechanical commutation process, 
which is analysed and measured as part of this research. 

The results of this research project have been written 
up in two major papers, the full texts of which are attached 
as the essence of this final report. These papers reflect 
the roost important ar.pects of the work. 

The first, entitled 

"A MODEL OF ANISOTROPIC ALNICO MAGNETS FOR FIELD 
COMPUTATION, ^ 

has been submitted to the I.E.E.E. for publication in the 
Transactions on Magnetics. It summarizes the text of the 
University of Southern California Ph.D. Thesis by Saad A. 
Al-Murshid, which has the same title. 

The second, entitled 

"COMMUTATION OF THE PERMANENT MAGNET AXIAL-FIELD 
D.C. MACHINE," 

has been submitted to the I.E.E. for publication in the 
proceedings part B. It summarizes the text of the University 
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of Southern California Ph.D. Thesis by Abdulhamid T. 
El-Kamodi, which has tha same title. 


A MODEL OF ANISOTROP(C ALNICO MAGNETS 
FOR FIELD COMPUTATION 

Peter Campbell and Saad A. Al>Murshld 
ABSTRACT 

The modelling of an anisotropic alnico magnet for the purpose of 
field computation involves assigning a value for the material's permeability 
in the transverse di: .,^tion. This is generally based upon the preferred 
direction properties, being all that are easily available. By analyzing the 
rotation of intrinsic magnetization due to the sslf<> demagnetizing field, it 
is shown that the common assiunptions relating the transverse to the 
preferred direction are not accurate. Transverse magnetization 
characteristics are needed, and these are given for Alnico 5, 5-7, and 8 
magnets, yielding appropriate permeability values. 
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INTRODUCTION 


In the calculation of magnetic fields in permanent magnet devices, 

modelling of the permanent magnet material itself usually relies upon data 

acquired from the major demagnetization curve. Most of the magnets used 

in electromechanical devices, such as rotating electrical machines, are 

anisotropic, and a curve measured in the preferred direction does not 

» 

represent magnetic properties in any other direction. Unless the 
permanent magnet is subjected to a uniform self- demagnetizing field, 
some information on transverse magnetic properties will be required to 
rind the complete field distribution. Of course, practical magnets do not 
approximate to the ellipsoidal shape required for a uniform internal field, 
and the boundary conditions, in addition to the magnet's anisotropy, will 
determine the extent of this field's non- uniformity. 

This paper will briefly review the approximations that are 
commonly made to account for the transverse properties of a permanent 
magnet, leading to a comparison with actual characteristics measured on 
some popular materials. By considering the analytical derivation of 
anisotropy, a simple technique is suggested by which the directional 
magnetic properties may be more accurately represented in a field 
computation. 

MAGNET FIELD EQUATIONS 

The demagnetization curve for a permanent magnet relates the 
components of B and H for the preferred direction, and these vectors are 
also related to the material's intrinsic magnetization M via: 


B « Uq (H + M) . 
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It is convenient to express the dependence of M upon H as 


and so 


M ■ ^ • 

1 * H + \XqMq . 


< 2 ) 

( 3 ) 


where = (1 + x) is an apparent relative permeability. 

For an anisotropic material, x and ja* are tensors, and may also be 
field dependent. However, the field quantities may be resolved into 
components in the preferred (p) and transverse (t) directions, so that 



and 

= 

This permeability tensor has been used to solve the field in regions 
containing soft>magnetic laminations [1], which are anisotropic but do not 
possess Mg. 

In permanent magnets, must be assumed to act only in the 
preferred direction, to allow and to fulfill transverse symmetry 
conditions in the material. Hence, 

M, . ] . 

However, this does not mean that the intrinsic magnetization is fixed in 

the anisotropic direction, rather that any rotation of will be accounted 

for bv X, and 

P t 

In an isolated circular magnet (Fig. 1), the self-demagnetizing 
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K will only lie In the p<>direction along the axis; elsewhere, there will 
also be a t- component of H, leading to a rotation of Mq that is governed 
by the material's anisotropy. The equations of the magnetization curves 
for these mutually perpendicular directions are deduced from (2): 


M a x, H + M 
P P P op 


Mt * 


(7) 

(8) 


The major curve (7) is a commonly available design characteristic, and is 
measured with a uniform field, H^, throughout the material. In contrast, 
little work has been done in the past to measure transverse properties. 

It has been assumed [2] that the permeability functions and 
are the same, and even [3] that is a constant, taking the value of at 
remanence. The permeabilities to be used will actually depend upon the 
operating conditions of the magnet. If these are static, the major curves 
are appropriate, and if the material was initially saturated in the 
p-direction, operation will be on the demangetization curve for the 
p-direction and on the initial curve for the t-direction; these conditions 
also satisfy the symmetry requirements setting = 0. These magnetic 
zation curves have been measured for Alnico 5-7 (Fig. 2) and Alnico 8 
(Fig. 3). Dynamic operation involves recoil, in the 2nd and 1st quadrants 
for p- and t-directions respectively, and a number of measured recoil 
loops are shoum in Figs. 4 and 5 for the same materials. A constant value 
for relative recoil permeability may be assumed throughout any quadrant, 
provided that only small transverse fields exist, and Table I gives a 
comparison of the values. It is apparent from the measured characteristics 
that, dependent upon the material type, Up and may be quite dissimilar. 
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Table I: Relative recoil permeabilities measured from 

Figs. 4 and 5. 



p-direction, 

t-direction. 


2nd Quadrant 

1st Qxiadrant 

Alnico 5-7 

1.82 

6.70 

Alnico 8 

1. 83 

2.89 


This use of the major magnetization curves implies that the p- and 
t- components are independent, which is in accord with (5). However, 
consider that in Fig. 1 there is a location at which only exists. This 
will rotate Mq to give a component as in (8), but (7) implies that 
remains with value This is a contradiction, and clearly the applica- 

tion of must reduce by some amount. Therefore, (5) does not 
provide a complete model of a permanent magnet, and we should include 

at least the off-diagonal term p' if not u* also. The more general 

pt tp 

permeability tensor is 



and we shall find values for all of these terms. 

PERMANENT MAGNET ANISOTROPY 

It may be shown that * p^p s 0 for magnets operating on a 

linear demagnetization curve, such as the ceramic ferrites and rare 

earth- cobalts. This linearity is derived from a constant , 

f P op 

meaning that Mq is truly fixed in the p* direction. It is the rotation of 
Mq that yields a non-linear curve, in materials such as the alnicos. The 
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general case of an Internal field H at angle 6 q rotating the magnetization 
by 0 is shown in Fig. 6, and the energy associated with this field 1st 

Ej « MqHMq cos (0q- 0) . (10) 

The torque that H applies to Mq is opposed by torques due to both 
shape and crystal anisotropy in alnicos. Considering that they have a 
cubic crystal lattice, the energies associated with these anisotropies are 
approximately given by [4]: 

Esh « K + K^j^sln^0 (ll) 

and 

E^r “ % 

where K, &nd are material anisotropy constants. The 

restoring torque on Mq is found from: 



where is the total of the energies in (10), (ll), and (12). Hence, 

T « -UqHMq sin (0Q-0) + sin 20 + sin 40 . (13) 

The usual technique for evaluating and is to rotate a thin 
disc of the material in a torque magnetometer, measurinf;: T and Bq. It is 
necessary to make 0 q ^ 0, and to eliminate the first term of (13), for which 
a uniform saturating field H is applied. 

The typical operating conditions of a magnet in a device do not 
involve a saturating magnetic field, and so measurements of T vs. Bg 
cannot be interpreted. However, Mg and 6 can be directly measured in 



a vibrating coil magnetometer [5], in which a moderate H is applied 
parallel to the plane of the disc. Two sets of coils, in line with and at 
right angles to the field, provide signals from which Mq is deduced, a» 
the sample is rotated through - 

A disc of Alnico 5 has been tested in this way, and we shall consider 
results obtained in a moderate field of 30 kA/m, which is somewhat less 
than the value required for saturation. Figure 7 shows {MqI vs. the angle 
of rotation of H» €q. and Fig. 8 shows the angular difference between H 
and Mq, 6q- € vs. 6 q. Shape anisotropy is predominant in Alnico 5, and 
Mq is much stronger in the p-direction than in the t- direction, with a 
preference for the axis of initial magnetization (6 qB 0). As the sample 
(\. e. , H) is rotated through 180°, stays closely aligned with H until 
the t-direction is approached, when Mq slips back almost to the preferred 
direction, and this is accompanied by the magnetization being nearly 
destroyed. However, as H approaches 180°, Mq reverses into alignment 
with H again, and the sample is pe.rtly remagnetized. 

These results are interesting as an illustration ox the rotation of 
Mq, but because of the large values of 6 that occur, they cannot be 
used to find the anisotropy constants. Even in a saturating field, Cq- 6 is 
still not small enough to counter the large value of H, to eliminate the 
first term of (13). Figures 9 and 10 show measurements at 211 kA/m, and 
now there is only a small variation in with 6^, the maximum 6 q> 6 being 
1, 2° at 0Q 45°. It will later be show'n that and Kj make the 
anisotropy energies small compared to the field energy here. 

Unfortunately, the use of a magnetometer for permanent x^iagnet 
materials does not yield the information needed to model their behavior, 
but we shall show that and Kj may be derived from the major 
denxagnetization curve. This is the same information from which u* 
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was inaccurately deduced for field computation [2, 3], and we shall suggest 
a more accurate model based upon the same simple test. As seen from 
Fig. 1. typical values for 6 q within a magnet are 150^ ^ ^0 ^ 210^. On the 
axis of symmetry, 6 q = 180^, and we shall find from the intrinsic demag- 
netization curve that « -53.7 kA/m, the field required to reverse M^. 
However, we also saw in Fig. 8 that a similar sudden change in the direction 

f 

of occurs with a lower field, 30 kA/m, directed at 6q ■ 117^. In 
general, a steady, predictable rotation of Mq will only occur within limits 
of { H| and 6 q, and if i 180^, the field at which a sudden change in 6 
occurs is less than 


ROTATION OF MACa^ETIZATlON 

An analytical model for a permanent magnet should include the 
effect that H has upon Mq. This may be found from (13), setting T s 0 
to find the equilibriiim position for Mq. Thus, 

M^HM^sin (0Q-6) * K^j^sln 26 + -g Kj sin 46 , (14) 


which may be rearranged as 




K 


sh 


[ sin (©Q-e) *] 

sin 26 I ’ 


K, 


K 


cos Z6 - 1 « 0 


sh 


(15) 


If the magnitude and direction (6 q) of H are known, and since (Mq| is being 
considered a constant, this may be solved give the rotation of which 
is 6. 

The critical field, at which the sudden change in 6 occurs, is 
found by setting [4]: 
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« 0 . 



Thus, from (13), 


H M.cos (6n~6) -f ZK . cos 20 f 2K, cos 46 « 0 
u c u u so 1 




cos 26 + 2Kj cos 46 

UjMgCO. (6g-e) 


_ 


( 16 ) 


For the major demagnetization curve, ^ci* ^0~ originally 

at 6 s 0, so 



2K . + 2K, 
sh 1 


07 ) 


It is useful to think of intrinsic coercivities due to the shape and crystal 

anisotropies separately, being and respectively. becomes 

H . if K, = 0. and H . if K . * 0. Hence, 

CIS 1 * cic sh * 


H . = H , + H , . 

Cl cis cic 


(18) 


and (16) may be rewritten: 

H . cos 26 H , COB 46 
u , CIS cic 

"c " cos (6 q-6) 

Thus, every angular rotation of may be considered to be unstable if 
the critical field, at 6^, corresponding to that 6. is experienced. 

It is more useful simply to know the value of corresponding to 
any 6q, for which 6 must be eliminated. The second expression needed 
is (14), a general relationship between H and Mq. which may similarly be 
rewritten as: 
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H 


( 20 ) 


2H . sin 26 4^ H , tin 46 
cts cic 

4 sin (6 q“6) 

The most familiar use of (19) is to draw the major demagnetization 
curve for 6 qS tt, 6 = 0. This theoretical characteristic is simply a square 
loop, with Mq reversing abruptly between 6=0 and 6 = n at± 
hysteresis loop for any other 6 q is plotted by calculating H from (20) and 
the component of M at 6^ from 

M = Mq cos (6q- 6) . (21) 

The continuous range of H is determined by (19)» Values for the constants 
are now needed. 

From a large block of the same Alnico 5, small 10 mm cubes were 
cut at 10^ increments (Fig. 11), and the intrinsic magnetization curves 
were measured (Fig. 12). Because of slight flaws in some of the cube 
samples, these curves do not exhibit an exactly steady progression. Mq 
is taken as the saturation or spontaneous magnetization, and ~ 1* 334 T. 

H . is measured from the major curve to be -53. 7 kA/m. According to (11) 

Cl 

and (^2), only crystal anisotropy gives a stable energy condition in the 
t-direction, so the 90^ demagnetization curve may be used to find 
H . = -30.6 kA/m, and hence H = -23, 1 kA/m. 

ClC CIS 

With these values, analytical intrinsic magnetization curves may 
be plotted for different 6 q as in Fig. 13. A study of hese curves Indicates 
any limitations in using (20) to model a real magnet. The most obvious 
discrepancy is the absence of remanence in the t-direction, which is 
because it was presumed that the spontaneous magnetization was in the 
p-direction. If the magnet were saturated with an field, the crystal 
anisotropy would then yield such a remanence. 
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It h*s b«en mentioned that reelietic conditions with a magnet will 
constrain H approximately to 150° ^ ^0 ^ 210°. In fact, there Is very 
good agreement between the measured and calculated remanence values in 
Figs. 12 and 13, from the 0° to the 60° curve, which encompasses 
120° < < 240°. Apparently the abrupt reversal of Mq does not occur 

in practice as predicted by (19), one reason being that the theory for shape 
anisotropy does not account for the interaction of neighboring magnetic 
elements in the non>magnetic material matrix. 

This method of modelling a permanent magnet from the rotation 
of Mq does not require testing of the material in a magnetometer, and 
will therefore be relatively simple to implement. It was mentioned that 
Kj and could not be evaluated in a torque magnetometer via (13), 
because the first term was not small. We now have these constants from 
Fig. 12, (17), and (18), which yield Kj » 20. 4kJ/m^ and « 15.4 kJ/m^. 
These are small compared to » 281 kj/m^ at H « 211 kA/m. 



PERMANENT MAGNET MODEI. 

We have already developed an analytical model that requires 
information from only the p- and t-direction magnetization curves, and it 
will now be shown how just these need be used to evaluate the terms of (9). 

Consider in Fig. 6 that H is applied at 6q ■ tt/2, for which (15) 
becomes: 


^o”t^o cos e 

*^sh 


K, 


K 


sh 


cos 26 


0 


( 22 ) 


The t- component of M^ is found from (21) to be 
Mj » Mq a in 6 


(23) 


10 > 


Eliminating 6 between these gives the transverse magnetisation curve, 

Mj vs. from vkrhich is found. However, a simpler function will result 
if it is assumed that only a small rotation of occurs, which is of course 
the case with anisotropic magnets in most practical applications. For 
small 6, 

K, 


1 


^0«t^0 

■ ^sh 


1 a 0 


(24) 


a M^e 


(25) 


Eliminating S, 


••• ^ ' Ch^ 


X. 2 


y “t 

Cl 


(26) 


(27) 


The rotation of will be reflected as a reduction of M^, such that 
-0 p 


“p = “o - 


••• M- = 


1 - 


iiti r 


Cl 


Expanding this, remembering that 1 Mq| = and approximating to the 

first term yields 


H 2 

M a M - -yM. (-TT^^ 
p op 2 OVH^. y 


(28) 


The permittivity lunctions may be written from (2) to include off-diagonal 
terms, as 

. x,H,+x^pHp 
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(29) 


Comparison between (27) and (29)* and between (28) and (30) yields: 


M, 


H 


ci 


’Sp 


pt 


= 0 


S 0 


*-c 


M 


2H : 

Cl 


Z J 


(31) 

(32) 

(33) 

(34) 


For the t>direction, the theory assumes that the material was 
spontaneously magnetized in the p- direction, and so alone will yield no 
M^, which leads to (32). The initial curve is modelled by (31)* which is the 
tangent to the 90^ curve in Fig. 13. This line coincides with the point at 
which saturation is achieved on the major (0°) characteristic. For the 
Alnico 5, the value of (31) is 19. 8, whereas the corresponding measured 
initial slope from Fig. 12 is 7. 3. The discrepancy is due to the measured 
value of = 53. 7 kA/m used in (31), and Fig. 12 shows that saturation 
actually occurs at a much higher field than this. There are two possible 
ways to achieve an accurate model of x^. 

If possible, the transverse magnetization characteristic should 

be measured, yielding directly, but this may be inconvenient since 

alnico magnets are difficult to grind, and their application most often does 

not require parallel, ground sides. Alternatively, there is a general rule 

in the magnetization industry [6] that the field required to saturate an 

alnico is three times its coercivity. Surprisingly, this works very well 

here, and the transverse curve is not needed. The coercivity is derived 

from the major curve in Fig. 12 as 50. 1 kA/m, and taking in (31) as 
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3x50. 1 kA/m yields s 7. 1, compared to the measured 7. 3. Either of 
these methods are a considerable Improvement over existing techniques 
[2, 3], and for comparison, the measured slope of the major intrinsic curve 
at remanence is only 0. 9. 

For the p- direction, the assumption of spontaneous magnetisation 
sets 0, whereas one might prefer to model the actual characteristic 
here as a function of Rotation of will cause the magnet to operate 
within the major hysteresis loop, as this process involves the loss of some 
magnet energy. Dependent upon the transverse field H^, will be 
reduced according to (34), but to take better account of the actual charac- 
teristics, this should be rewritten as 


s - ( — H 

*pt V 2M. y ”t 


(35) 


Using the values found for and Mq for Alnico 5, a field of = 30 kA/m 


will cause only a 0. 3% reduction in M^. This is extremely small, and 
leads to the conclusion that although x^^ is not actually zero, it may be 
assumed to be so. 

The apparent relative permeabilities equivalent to these x- functions 
are written by combining (29) and (30) with (1). The analytical form of (9) 
then becomes 


\L' * 


1 


2M.. 


1 + 


(36) 


However, in the light of the experiments described, it is preferable to 
use the measured p- and t- characteristics directly, in which case (9) 
reduces to its approximate form, (5). 
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Different investigators will have their own preferred methods for 
modelling the major demagnetisation curve, with as a function of if 
necessary. Generally, though, information on transverse properties is 
not available, and we conclude with some further data obtained by grinding 
cubes of Alnico 5-7 and Alnico 8. Their magnetization curves were shown 
in Figs. 2 and 3 respectively. We have stated that symmetry requirements 
dictate the initial t- characteristic be used for static field applications, and 


this is approximately linear in small fields. The constant values of that 
should be used in (5) are given in Table II. 


Table II: Relative Initial Permeabilities Measured in the 

Transverse Direction. 


t 



Material 


: ( 

Alnico 5 

8. 3 


Alnico 5-7 

9.25 


Alnico 8 

1.9 




CONCLUSIONS 

Whereas it is common practice to set off-diagonal terms in the 
permeability tensor to zero, experimental data and analysis have been 
presented that justify this assumption. Moreover, important information 
on the transverse characteristics of the three most common permanent 
magnets will help to improve field computation. The values of for 
Alnico 5 and 5-7 are particularly high compared to as was the case 
for recoil permeability, and this would tend to enhance transverse flux 



I . 

i 



In a magnet that directly experiences armature reaction in an electrical 
machine. If an exact model of a permanent magnet is required, then the 
rotation of is found using (15) or (20). 
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Fig. 10: Angle by which H leads Mq » with rotation 

by 6^ of H « 211 KA/m, for Alnico 5. 










Fig. 13: Intrinsic magnetization curves calculated at various 

angles 0^ to the preferred direction in Alnico 5, 
using y M = 1.334 T and H . =* — 53.7_KA/m. 
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ABSTRACT 


The power losses In a permanent magnet d.c. disc-armature motor are 
analysed, and the predominant component Is found to be due to short-circuit 
currents In colls undergoing conmutatlon. Measurement of these currents 
shows their magnitudes to be greater than predicted from an equivalent 
circuit of the machine, and this suggests that the brush contacts may not 
be modeled by resistance alone. Inductances used In the equivalent circuit 
are very small, and their angular variations are found to depend upon the 
properties of the anisotropic alnico magnets; the transverse recoil permea- 
bility of Alnico 5-7 Is about four times that In the preferred direction, 
giving some sallency. 
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e.m.f. In one winding path. 

e.m.f. In a single coil. 

short-circuit current in a coil. 

mesh current in neighboring winding paths. 

contact resistance between one brush and coimutator. 

contact resistance between one brush and one coninutator segment. 

coil resistance. 

time at onset and end of conmutation. 

vector magnetic potential. 

flux density. 

armature current. 

winding current density. 

coll self-inductance. 

total self-inductance of two neighboring winding paths. 

total mutual inductance between short-circuit coil and two 
neighboring paths. 

mutual Inductance between short-circuit coil and first two 
neighboring colls. 

commutation power loss. 

constant term of power loss. 

resistance added In series with coils. 

total machine resistance. 

angular displacement between adjacent coils. 

relative permeability. 

angular position of short-circuit coll. 

armature angular velocity. 

total flux per pole. 



1. Introduction 


The permanent magnet d.c. disc-armature motor has a well-known 
topology,^ a heteropolar axial-field magnet structure being adjacent 
to an Ironless armature, the winding being encapsulated In epoxy resin. 

The characteristically long air-gap gives, this machine a low armature 
Inductance and a relatively high field about the magnetic neutrals. Both 
of these are expected to result In high power losses associated with the 
mechanical conmutatlon process, and experiments on an early disc-armature 
motor have supported this view. However, data acquired on that machine 
was limited because It was not possible to modify the armature for direct 
reading of currents. 

The addition of resistors In series with one end of each coll and the 

3 

comnutator bars provides direct measurement of coll current via slip rings. 

A variation on this technique Is used with a specially constructed disc- 
armature n»tor, to verify computed short-circuit currents In a coll under- 
going commutation, and thence to find the contribution that these make to 
the power losses In the machine. The experimental motor has eight poles, 
and a 72 coll lap-connected winding. After conventional steady-state load 
tests were performed, the armature was modified by the addition of the 
series resistors. Although only one resistor Is connected to the slip 
rings for coll current measurement, a balancing resistor must usually be 
placed In all other colls. ^ However, the total addition to a winding path 
(9 colls) may adversely affect the performance of the machine. 

The technique employed In this Investigation Involves the addition of 


( 1 ) 


just one resistor to each parallel path, the value of 0.12n having little 
effect on performance and yet providing an adequate measurement signal. 

The armature of the experimental machine has an outer active diameter of 
221imi (Fig. 1), and space Is left between winding and coonutator to accom- 
modate the resistors and silver slip rings. Rather than having eight 
resistances, one added to every 9th coll, there being 9 colls/pole, there 
are nine resistances, one added to every 8th coll. The machine was con- 
structed with brush arc and commutator segment arc equal, so that just one 
coll Is connutated by each brush at any time. Therefore, it any time each 
parallel path has one added resistance, and one coll with a resistor Is 
being connutated (Fig. 2). Balanced parallel winding paths are maintained 
while short-circuit coll current Is recorded. 

2. Equivalent circuit of the dlsc-annature motor . 

An analytical technique Is described for the evaluation of coll 
current and the associated power loss of coiimutatlon, for the axial-field 
permanent magnet d.c. machine. This will be verified by test data from the 
experimental motor. 

The appropriate equivalent circuit for a coll undergoing connutatlon 
Is shown In Fig. 3. This coll has resistance r^, self-inductance L, an;2 
encounters an Induced e.m.f. e^. The particular coll under study will also 
have the added series resistance R, and the short-circuit current will flow 
through one resistive contact r^, to the brush, returning through another 
contact r^^^. The circuit Includes the neighboring amature paths, each of which 
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shares mutual Inductance VI with the conimitatlnq coll. Since this Is a 
lap winding with an exact number If colls per pole, each neighboring path 
has approximately the same e.m.f. e^. self-inductance and resistance 
8r^ 4 R (g colls/pole less one being conniutated). No field winding appears 
In tne equivalent circuit, but If the permanent magnet material Is an 
alnico type, there will be a salient-pole structure which will cause the 
Inductances to be dependent upon the armature's angular position. Because 
of synmetry between the neighboring paths, each parameter In one will be 
equal to that In the other. 

By suitable choice for the mesh current 1^ In Fig. 3, the terminal 
voltage V Is eliminated and the e.m.f. s e^ cancel. The two mesh equations 
In 1 and 1^ are 



These may be solved for the short-circuit current 1, provided that all of 
the resistances and Inductances are known. 

2.1. Resi Staines. 

The resistance of a winding path Is simply measured by a static drop- 
test at rated current, from t«h1ch r^ Is deduced as 0.042n. Brush contact 
resistance, however, will vary both with angular velocity and angular position 
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of the armature. The study of brush contact phenomena Is alone a most 
Involved subject, and a simplified representation was employed for modeling 
this machine. It was assumed that the contact resistance between a brush 
and a conmutator segment is inversely proportional to the contact area, and 
that r^^ and r^^ (^^9* with angular position ut according to 


bl 


* r. “^2 * *^^1 ; r, 


(i)t > ii>t. 


b2 


1 


_ Ut- - (i)t< 

r^, 1 

° wt2 - ii>t 


(3) 


where t^ < t< t^ &nd t 2 give the times at the onset and end of conmutation 
for the coil under study. The parallel combination of r^^ and T |^2 * 

which is independent of t; this is the effective resistance seen through a 
brush, which may be deduced from load tests on the motor. This procedure is 
described in sec. 3, from which the variations in r^<j and r ^2 found with 
the aid of Eq. (3). 


2.2 Coil e.m.f. 

To obtain the greatest power output from the experimental machine 
with the given armature, the highest possible air-gap flux density should be 
obtained. This is achieved with high remanence Alnico 5-7 permanent magnets, 
assembled in two sets, one on either side of the armature. This arrangement 
is preferable to having a single magnet set to one side and a steel ring on 
the other, as described in Ref. 2, as it reduces the leakage flux. 


Because the ironless motor of the disc-armature machine ensures that 
negligible armature reaction occurs, a slow measurement of flux linking on 
actual coil may be used to determine the e.m.f., e^. The flux variation was 
found to be almost exactly sinusoidal, which is typical of this type of 


( 4 ) 



motor with & large air-gap. The experimental machine has 8 poles, the peak 
flux linking a coll Is a « 10.55 mWb, and so the coll e.m.f. Is given by: 

e^ • 4«i»a.cosamt 

« .O42u.cos4mt. (4) 

The annature angular velocity Is m (rad/sec), and connutatlon occurs 
about the magnetic neutrals, corresponding to ut « x/8, 3ir/8, 5ir/8, etc. 

The salient pole structure and the much smaller air-gap length of conventional 
small d.c. machines tend to supress the coll e.m.f. somewhat during commuta- 
tion, which is In contrast to the rapid Increase In e.m.f. about the neutrals 
Indicated In Eq (4). Acting In conjunction with low inductance values, this 
Is the cause of significant short-circuit power losses In the disc- armature 
motor. 

2.3 Inductances, considering Permanent Magnet Properties . 

Calculation of the self and mutual inductances of an armature coll 
Involves evaluation of the field resulting from rated current In the coll. 

The permanent magnets are Immediately adjacent to the armature; If these 
were rare earth-cobalt or ferrite magnets, they would appear to the coll as 
air-gap, but alnicos have non-linear demagnetization characteristics. The 
relative permeability, u. Is variable along the B-H curve of an alnico and 
the recoil permeability Is also different from unity. The consequent salient- 
pole structure gives all the Inductances their angular dependence. 


Many alnico magnets are also anisotropic, as seen by the major 6-H 


loops of Alnico 5-7 in the preferred and transverse directions (Fig. 4). 
These curves were measured on the actual material used for the experi- 
mental machine, one magnet being ground to a cube for testing In a 
permeameter. The conventional method of modeling such magnets that are 
subjected to an external field, as by armature coil, is just to use the 
recoil permeability from the second quadrant of the major B-H loop. Some 
recoil loops for the Alnico 5-7 magnet are plotted in Fig. 5, for which a 
constant recoil permeability of 1.82, the mean value, may be taken as a 
good approximation. 

However, it is well-known that field within a non-linear permanent 

5 

magnet may not be constrained to the preferred direction, and so transverse 
components must be considered. In this case, the first quadrant of the 
magnetization loop should be usectwith recoil coming from the virgin curve. 
This is to satisfy the necessary conditions of synmetry about some plane 
within the magnet where no transverse field exists. Figure 6 gives the 
measured recoil loops for the Alnico 5-7 magnet, and since only small 
transverse fields will generally occur, a constant value of 6.70 may be 
taken for the recoil permeability. 

The recoil permeability used to model Alnico 5-7 is therefore almost 
four times greater in the transverse direction than in the preferred 
direction. Remembering that alnico permanent magnets are generally many 
times longer than the air-gap, and that the steel backing rings are well 
removed from the region, it is possible that the self-inductance of a coil 
may be greater with its axis over the magnetic neutral than with it over 
the magnet center-line. This is contrary to the situation in a conventional 
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salient pole machine. 


Since the active conductors In the axial-field machine lie radially, 
an approximate field distribution In the air-gap region may be found 
from a numerical solution of Poisson's Equation In 2-d1mens1ons, at a 
fixed radius. The appropriate form of this will take account of the 
directional dependence of the permeability u In the magnets, so that the 
flux density 6 is calculated via vector potential A from the current 
density distribution using 

V X B ■ A (5) 

and ^ ^ 


Figure 7 shows a typical field distribution calculated at the average 
active radius for rated current In a coll, whose axis Is midway between 
the magnetic neutral and the magnet center-line. The effect of having 
the higher transverse permeability Is to shunt magnets In the vicinities 
of the conductors. Here we are evaluating Inductance; however. It Is 
fortunate that armature reaction Is small In this machine, because other- 
wise this effect would seriously complicate any attempt to calculate the 
resultant (permanent magnet plus armature reaction) field distribution. 

Account was also taken of end-winding effects, but since these are 
well removed from any magnetic materials, their contribution to Inductances 
was found to be small. The angular variations of Inductances (Fig. 8), 
determined from calculated field distributions, do verify that the 
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principal Inductances are greatest with the coll axis over the magnetic 
neutral regions. The se1f>1nductance L for this 8 pole machine may be 
written as the series 



L • Lq + cosSut + L 2 cosl6wt + cosSZut, (7) 

where Lq * 8.06yH; « 1.76yH; L 2 “ -.46yH; * .16yH. 

Mutual Inductance with a neighboring coll depends both upon 6 ■ lat, 
the position of the source coll, and the relative position between the 
pair of colls. There Is a strong coupling of the sinusoidal variation 
with the Immediate neighbor, M2, and also with the coll that Is exactly 
one pole pitch away, M10. 

Measurements of Inductance were obtained by constructing an armature. 
Identical to that In Fig. 1, but without a comnutator and hence without 
Interconnections between colls. Assembled In the stator of the experi- 
mental machine, with one coil excited at rated current, the actual varia- 
tions of Inductance with angular position were found directly (Fig. 9). 
These all have extremely low values, characteristic of the Ironless rotor 
permanent magnet machine, and are generally less than the calculated 
figures. The discrepancy Is due to having taken a two-dimensional field 
approximation, neglecting the effect of the magnet boundaries at the inner 
and outer active radii. A three-dimensional field computation should 
yield more accurate results.® 

For the parameters to be used In Eqs. (1) and (2), It Is assumed 
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that L, M2 and M3 (the first two neighbors), have purely sinusoidal 
variations with ut. higher hannonics being Ignored; mutual Inductance 
with all other colls Is taken as constant. From the results of Fig. 9, 


L • Lq 4 cosSuit 

(8) 

M2 “ M2q 4 M2^ cos8((at 4 ijy) 

(9) 

M3 « M3q 4 M3^ cos8(ii)t 4 -y) 

(10) 


where Lq • 7.785yH; « .ZESpH; 

M2q * 3.64pK; M2^ « .SlyH; 

M3q » 2.61yH; M3^ * .24yH. 

If two adjacent colls are displaced by an angle y, being 5^ In this 
72 coll winding, the peak mutual Inductance between them will occur at 
%Y from the position corresponding to the peak self>1nductance (Eq. 9). 
Similarly, with M3 the variation Is shifted by 2 x Vr In Eq. (10). 

The total self>1nductance, of the sixteen colls In the two 
neighboring armature paths In Fig. 3 Is deduced from Eq. (8) to be 

8 

* 16Lq + 2L^I cos8((Dt ♦ ny) (11) 

n*l 

Similarly, the total mutual Inductance linking these paths to the 
short'd rcul ted coll Is found using Eqs. (9) and (10) to be 

9 

M • 2l MNft ♦ 2M2, cos8(ut ♦ *fy) + 2M3, cos8(wt 4 y) (12) 

■ic9 V I • I 
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The first term Is evaluated from Fig. 9 as 47.3yH. 

An important conclusion of this analysis is that, when an electrical 
machine topology involves a1nico*type permanent magnets adjacent to the 
armature, consideration must be given to the magnetic characteristics 
perpendicular to the preferred (anisotropic) direction. Calculation and 
measurement of inductances has demonstrated the importance of field 
traversing the magnet faces. In Sec. 4, the calculated short-circuit 
current during Conmutation is compared favorably with measured variations. 

3. Determination of Motor Parameters and Performance . 

0 

It was stated in Sec. 2.2 that the brush contact resistance will vary 
with angular position according to Eq. (3), averaging to a value of 
at each brush; it will also be dependent upon the motor's speed. It will 
therefore be necessary to perform a load test on the motor at a constant 
speed, to determine the resistance that is associated with that speed. 

The experimental disc-armature motor is shown in Fig. 10, the armature 
of Fig. 1 being located between the two sets of Alnico 5-7 magnets. In 
addition to the added resistors (R) in the winding, there is a 0.20 
resistance (mounted above the machine in Fig. 10) in series with each 
brush lead for measurement of the brush current. The windings on the 
permanent magnets are for the purpose of initial magnetization only. 

• 

Rated voltage and speed for this machine are 60V and approximately 

2400 rev/min. A load test at this constant speed yields the relationship 

2 

between the total power losses and the square of armature current (I ) 
given in Fig. 11. Except at very low current, the losses may clearly be 
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represented by a constant term, • 192W, and a constant resistance, 

« .309R. P|^ Is conventionally treated as a mechanical loss torque, 
and accounts for all the electrical losses at this speed. 
An expression for the power losses involving output torque, T, Is: 

VI - mT • I^R^ + (13) 

Considering the values for the Individual resistances In this eight 
pole lap-wound machine, given here and In Sec. 2.1, the total resistance 
of ■ .309n breaks down Into the following components: 


(1) armature winding resistance. Including additional R In 


each path: .057Q 

(II) total brush contact resistance: .1520 

(III) external resistances .1000 

Total : .3o9& 


At each brush contact, we deduced a value for r^^ » .3040, which Is to be 
used In Eq. (3) to solve for the short-circuit current from Eqs. (1) and 
( 2 ). 


The complete performance characteristics on this motor at 60V are 
given In Fig. 12, and over the measured current range, speed Is seen to 
vary between rated speed and IIS less than this. Extrapolation of the 
speed curve to a very high stall current Is characteristic of permanent 
magnet d.c. motors. With such a small variation In speed over this 
current range, we would expect good correlation between the power losses 
here and those at 2400 rev/mln In Fig. 11. 
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This is so, and the value for derived from Fig. 11 may be used wlth^ 
confidence to evaluate motor performance at GOV for 0 < 1 ^25. 

With the added resistors In the armature winding and the brush leads, 
the current was held below Its normal rated value of 40 A for these tests. 
Heating losses In these components cause the efficiency in Fig. 12 to 
peak at about SOX. This figure Is higher without these added resistors, 
and performance characteristics obtained on the machine prior to modifi- 
cation of the armature (Fig. 13) show a peak efficiency of 851. In this 
case, the armature winding resistance was .0420, and the brush contact 
resistance of .1520 represented 78X of the total resistance. This high 
proportion Is Indicative of a problem occurring at the brush contacts, 
and this will be shown to be the source of the dominant power loss In 
this machine. 

4. Conroutatlon Power Losses . 

The objective Is to compute the power loss associated with short- 
circuit currents during cormtutatlon, and to determine whether they appear 
In Fig. 11 as part of the constant term, or the variable 1 loss. 

4.1 Short-Circuit Currents during Commutation . 

Equations (1) and (2) are solved for the mesh current 1^ (Fig. 3), 
;«h1ch will be zero at the onset of commutation for the short-circuited 
coll, and will only vary through the mutual Inductance with 1. This Is 
due to the path chosen for 1^, and hides the fact that direct terminal 
current flows In both of these neighboring paths, through the brush 


resistances and r^. Clearly, this was taken Into account when 
deriving values for r^^^ and r ^2 constant speed test, and these 

resistances should now be used with the appropriate currents. This is 
achieved by setting the initial and final values of 1, at the onset and 
end of conmitation, equal to the normal path current. 

As an example, a terminal current of 24. 4A will divide Into the 8 
parallel paths, giving initial and final values for 1 of -3.0SA and ‘*^3.05A. 
At the onset of comnutation for the coil In Fig. 3, it has Just emerged 
from the right-hand neighboring path, and the initial current of -3.05A 
accounts for the losses in the r^^^ and r^^^ contacts in that path. The 
power losses in the contacts in the left-hand path are equal to this, 
being accounted for by another coil entering comnutation at the next brush. 
There are, in fact, eight coils in the machine in the same condition cor- 
responding to the eight parallel paths. 

The calculated variation of current in a coil undergoing commutation 
is shown in Fig. 14, for operation at 60V, 24. 4A, and a speed of 2171 
rev/min derived from Fig. 12. Because the armature coils have very low 
Inductance values, comnutation is achieved very rapidly. Mutual coupling 
will cause a spike to appear on the current of the neighboring coils, and 
conversely, spikes appear on Fig. 14 at the instants when neighboring 
coils are comnutated, these spikes should accumulate and appear on the 
terminal current waveform, and Fig. 15 shows the measured current 
entering one positive brush. The variations in the spikes are due to 
non-unifonnities In the armature winding and the magnets. 
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The measured variation of coll current Is shown In Fig. 16, for the 
same operating conditions. These oscillograms confirm the predicted rapid 
reversal of current at the onset of commutation, but show It rising to a 
higher peak than calculated during the first half of the conmitatlon period. 
One possible cause of such a discrepancy would be If the brushes could tilt, 
connutatlon would be advanced and the peak current would Increase. However, 
solution of Eqs. (1) and (2) showed that this effect would not yield the 
desired peak current. 

It Is postulated that the brush contact Is not simply resistive, but Is 
also capacitive. In that case, the path for the short-circuit current, 1, 

In Fig. 3 forms a series resonant circuit, which would requlrr a capacitance 
of about 480 yF at a resonant frequency corresponding to 2171 rev/ml n. If 
such a capacitance existed In a conventional machine, the normally much 
higher coll Inductance would reduce the resonant speed to a very much lower 
value, and an effect such as this would not be detected. It Is not expected 
that such an effect will have a simple relationship with terminal current, 
but It should still be evident for motor operation on no load. Figure 17 
shows the coll cur**ent at 60V, 3.2A (no load), and Indicates that the high 
current peaks do still exist, but are variable and somewhat unpredictable; 
this Is typical behavior of brush contact phenomena. To support the 
theory that capacitance exists, close examination of the coll current on 
load (Fig. 16 (b)) shows that ringing occurs after the reversal of current 
at cofnnutatlon. 

4.2 Fower losses. 


From the calculated variation of 1 In Fig. 14, the power loss In the 


coll can be evaluated for the entire comnutatlon period. A similar 
variation, without the added resistor R In the coll, will yield the power 
In the other seven colls being conmutated. The effect of terminal current 
has been accounted for In the Initial and final values that were set for 
1, and so these eight colls represent the total comnutatlon power loss. 

Some calculated values are given In Table I, together with the power loss 
In the other fixed resistances In this machine. 

The total calculated electrical power loss at 2400 rev/mln Is sub* 
tracted from the measured curve In Fig. 11. The mechanical power loss 
at this spebd was measured by demagnetising the stator of the experimental 
machine, and by driving It with another Iron-less axial -field motor. There 
remains a further component of the total losses In Fig. 11. which Is also 
associated with the comnutatlon process. If this component Is added to 
the proportion of the calculated electrical loss already attributed to 
commutation, then the comnutatlon power loss Is represented by: 

• 1^5 4 .152 (W). (14) 

comn 

Comparing this to the total power losses given In Eq. (13), commutation 
accounts for 76X of the constant component, P^, 49X of the current depen- 
dent component In the experimental machine, which would be 78% of r*R^ In 
the machine prior to modification of the armature. 

That the comnutatlon power loss Is significantly* Independent of 
armature current. Is supported by the oscillograms of short-circuit 
current, which show that the peaks at the onset of .comnutatlon on load 
(Fig. 16} are still present on no load (Fig 17). 

(15) 
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Table I: Calculated I R power losses for 60V supply, using r. measured at 
2400 rev/ml n. 


Terminal 
current, 
I (A) 

Total 

connutatlon 
power 
loss (W) 

I^R loss In 
winding and all 
added resistors 
(W) 

Total I^R 
power loss 
In motor 
(W) 

3.2 

8.1 

1.6 

9.7 

18.0 

87.7 

50.9 

138.6 

24.4 

154.0 

93.5 

247.5 




Armdture re&ctlon, not present In this nachlne, would cause this power 
loss to be more heavily dependent upon terminal current. Although the 
total power losses may be represented by Eq. (13) » from which a steady 
value for r^^ was deduced, this resistance does not have to account for 
the entire brush contact loss. That Is, the brush contact does not have 
to be simply the constant r^ as seen through the terminals. The evidence 
Is that the contact should be modeled also by a re:>ct1ve tenn, chosen 
such that the circulating current matches the oscillograms (Figs 16, 17). 

S. Conclusions. 

« 

Despite having a large proportion of Its losses attributed to unwanted 
circulating currents In colls during coomutatlon, this experimental disc- 
armature motor exhibits a peak efficiency of 85% at rated current, and 

maintains above 80% efficiency over a wide current range from about 15A 
to several times rated current (Fig. 13). This coninutatlon power loss Is 
evident even without load, and must be considered to be the majcr problem 
In obtaining optimum performance from this type of low Inductance machine. 
The Ironless armature gives benefits such as elimination of rotational 
Iron losses and armature reaction, but these must be weighed against the 
effects of having significant field In the magnetic neutral regions, 
creating the e.m.f. that drives the short-circuit currents. Under current 
Investigation Is the Inclusion of magnetic wedges between coll sides In 
the armature, to enhance and constrain the field under the poles, these 
wedges being made from a low loss Iron powder composite.^ 

Since coniDutatlon power loss Is the critical problem In this machine, 
analysis should b;; performed using an equivalent circuit as described here. 

(16) 


However, it Is apparent from the oscillograms that Fig. 3 may not be 
complete, because the short-circuit currents and power losses were found 
to be predominantly independent of terminal current. Therefore, whatever 
impedance is used for the brush contact should be derived from measurements 
of the circulating current, rather than the terminal current. Further 
and more detailed studies using this technique to measure i will be re- 
quired to determine that impedance. 

The very low values for the inductances in this machine cause 
comnutation to occur very rapidly, leaving a coil short-circuited for a 
much longer' period than is necessary simply to reverse the current. The 
benefit of having higher inductances appears to be that they would eliminate 
any resonance that is causing i to rise to such high measured values. 

This paper has described the calculation of inductances in this permanent 
magnet machine, and so the effect upon performance of adding any magnetic 
material to the armature could be computed. In this case, the properties 
of the permanent magnets adjacent to the air>gap will be important, and 
for the anisotropic Alnico 5*7 described here, the recoil permeabilities 
in the preferred and transverse directions are qu»te different. Conmu- 
tation occurs with the coil and pole axes coincident, the position of 
minimum coil self-inductance in this type of machine. Current research 
is aimed at broadening the investigation to other important types of 
permanent magnet material, in particular Alnico 8, to compare the effect 
of different field distributions upon the short-circuit currents and power 
losses. 
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Fig. 5: Recoil loops in the second quadrant measured for Alnico 5-7 in 

the preferred direction. 
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Fig. 6: Recoil loops in the first quadrant measured for Alnico 5-7 in the 

transverse direction. 
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Fig. 8: Calculated variations of the self-and mutual Inductances with the 

source coll axis 0® from the magnetic neutral. 





Fig. 10: The experinentel di sc-arnature moior, with one set of Alnico b ’ 
maoncls either side of the arr.alure, showing resistors connected in the 
brush leads. 
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Fig. 11: Total power loss measured at constant speed • 2400 rev/mln for 
the experimental machine, after modification with added resistors, showing 
components of ti^ losses. 
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Fig. 12: Perfom.ancc curves at 60V for the experimenUl machine, after 

modification with added resistors. 
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ig. 15: Current entering one brush, for the machine operating at 60V, 

4.4A, 2171 rev/min. 
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rig, 16: Measured variation of the coil current, i, for the motor operating 

at BOV, 2A.4A, 2171 rev/min. 
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